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Abstract —  Results  from  design,  synthesis,  and  analysis  of 
optimal  mutually  dispersive  symbols  are  presented  as  an 
improvement  over  existing  symbol  sets  employed  in  Non-linear 
Ambiguity  Suppression  (NLAS).  A  recently  proposed  theorem 
formulates  the  existence  of  symbol  families  having  optimal 
mutual  dispersion,  a  highly  desirable  property  for  NLAS 
applications.  Results  from  theorem  analysis  are  presented  and 
compared  to  other  suitable  NLAS  symbol  sets,  showing 
significant  improvement  in  mutual  dispersion  characteristics. 
NLAS  ambiguity  suppression  effectiveness  is  demonstrated 
using  a  set  of  optimal  mutually  dispersive  symbols. 

I.  Introduction 

Radars  employing  pulsed  waveforms  are  inherently 
ambiguous  in  range  and  Doppler.  In  1962  Palermo  [1]  used 
two  conjugate  linear  frequency  modulated  (LFM)  pulses  to 
demonstrate  a  Non-linear  Ambiguity  Suppression  (NLAS) 
signal  processing  technique  for  reducing  ambiguous  energy 
in  processed  radar  returns.  The  use  of  conjugate  LFM 
pulses  for  diverse  pulse  coding  does  not  extend  to  larger 
symbol  families  and  thus  has  severe  limitations  for  M- 
channel  (M  >  2)  NLAS  applications.  Desirable  NLAS 
symbol  sets  posses  1 )  large  partial  period  autocorrelation 
peaks  with  low  integrated  sidelobe  levels  and  2)  maximum 
signal  dispersion  when  cross-correlating  pulse  codes  within 
the  family  [2],  Symbol  sets  comprised  of  pseudo-random 
discrete  codes,  such  as  Gold  codes,  exhibit  mutual 
dispersion  properties  [3]  but  are  not  optimum  for  NLAS.  A 
root  mean  square  (rms)  time  duration  metric  of  correlation 
functions  is  introduced  to  formulate  a  process  for  obtaining 
optimal  mutually  dispersive  NLAS  symbols  for  arbitrary 
code  family  sizes  [4]. 

II.  Symbol  Design 

The  discovery  of  symbol  families  having  desirable 
characteristics  for  NLAS  applications  is  by  no  means  a 
recent  achievement.  Code  Division  Multiple  Access 


(CDMA)  communication  schemes  use  discrete  codes 
exhibiting  desirable  properties  for  NLAS  [3].  The  search  for 
optimal  symbol  sets  has  historically  focused  on  discrete 
codes;  an  avenue  of  research  not  yet  providing  an  optimal 
solution.  Using  rms  time  duration  of  correlation  functions 
as  the  metric  for  optimality,  the  search  for  mutually 
dispersive  codes  begins  by  considering  a  pulse  coded  radar 
signal  given  by: 

M  - 1 

x(0=  Ifk(t-kT)  (1) 

k  =  0 

where  T  is  the  Pulse  Repetition  Interval  (PRI)  and  fk(t)  has  a 
Fourier  transform  of  the  form: 

Fk(co)  =  ^A(co)ej0k(co)  A(co)  ,&k(co)e3l  (2) 

The  processed  pulse  return  is  represented  as  follows: 

*  4 

Fi  ( co)Hk  ( co  )  =  Fl  ( co)Fk  ( co  )  =  Plk(co) 

=  A(co)ejf0l(a,)-0k(co)] 

When  k  =  /,  the  processed  output  represents  a  matched  fdter 
response  whose  inverse  Fourier  transform  is  the  signal 
autocorrelation  function  (ACF),  pkk(t)  given  by: 

3 -I 

pkk(°>)  =>  Pkk(0 

The  NLAS  process  requires  the  ACF  to  be  as  compressed 
(impulse-like)  as  possible  and  the  crosscorrelation  function 
(CCF)  (mismatched  filter  response)  to  be  as  dispersed  (flat) 
as  possible.  The  rms  time  duration  of  correlation  functions, 
Oik  as  defined  in  (5a),  is  used  to  quantify  correlative 
dispersion  [5].  This  definition  of  <J/k  assumes  P/k(t)  is 
normalized  to  unity  energy  with  zero  mean;  the  dot  operator 
(•)  in  (5a)  and  (5b)  represents  differentiation. 
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For  k  =  1  in  (5b),  crtt  represents  the  ACF  rms  time 
duration  and  is  purely  a  function  of  envelope  /)  t fty  (second 
term  of  second  equation  identically  cancels).  To  achieve  a 
maximally  compressed  ACF  response,  (5b)  is  made  as  small 
as  possible;  this  provides  the  impetus  for  designing  optimal 
envelope  A(co)  such  that  the  ACF  rms  time  duration  is 
minimized.  Similarly,  to  achieve  a  large  dispersed  response, 
the  CCF  rms  time  duration,  k  ^  I  in  (5b),  is  made  as  large  as 
possible;  the  aik  dispersion  has  contributions  due  to  both  the 
envelope  and  the  phase  functions.  With  the  envelope 
minimization  constraint  in  place,  the  phase  functions  are 
designed  to  optimize  the  CCF  rms  duration. 

Optimality  in  compression  and  mutual  dispersion  is 
achieved  using  the  following  process  [4]: 


(9 


1  2  1  k  =  l 
—  \(pk(co)(pi(co)A-(co)dco  =  8kt  =  \ 

2k  [O  k*l  ) 

4)  Describing  M phase -rate  functions  as: 

®k((0)  =  jGDck<P(a>) 

where  GD  is  a  dispersive  gain  factor,  (dot  represents 
differentiation)  yields  phase  functions  of  the  form: 


5) 


^k(a)  =  yfGo  \(p(co)dco  (10) 

CO 

Resultant  time  and  frequency  domain  expressions  for 
mutually  dispersive  codes  is  expressed  as: 


fk(t)  =  jA(^)ej[(0t  +  (S>k(a)]dco 

2n  -no/ 2 


(11) 


Fk(co)  =  V  A(o> )  ej°k(03) 

It  can  be  shown  that  symbols  designed  with  this  process 
yield  theoretical  ACF  and  CCF  rms  time  durations  of: 


1)  For  optimal  compression,  the  envelope  taper  is  derived 
by  solving  solution  the  following  constrained 
optimization  problem: 


J  •  2 

Minimize  —  \A  (co)d(O 

1  °°  2 

Such  that  —  \A  ((o)dco  =  1 

2n 


(6) 


Using  calculus  of  variations,  a  cosine  taper  envelope  of 
the  following  form  may  be  derived: 
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Furthermore,  these  theoretical  expectations  are  independent 
of  basis  selection  provided  the  aforementioned  process 
requirements  are  satisfied. 


III.  Symbol  Synthesis 


A(  co)  = 


\co  <- 


Qn 


(7) 


2)  For  M  symbols,  use  any  set  of  M  maximally  equidistant 
unit-vectors  in  M-Dimensional  space  (hermits), 
C=[c0  Ci  Ci  ...Cam].  One  solution  to  this  hermit 
problem  is  defined  as  follows  [4]: 


C  = 


1 

Jm2 -m 


M  —  1 

-1 

-1 

-1 

M  —1 

••  -1 

-1 

-1  ■■ 

-1 

-1 

-1 

...  M- 

IMxM 


(8) 


3)  Choose  any  set  of  odd  basis  functions  given  by 
cp(co)  =  [cp0(co) ,  cpdco) ,  cp2(co) ,  ...  ,  cpM_i(co)f 
while  satisfying  windowed  orthonormality  given  by: 


A  closed-form  expression  for  symbol  fk(t)  is  generally 
unobtainable.  Consequently,  a  rectangular-rule  approxi¬ 
mation  to  the  inverse  Fourier  transform  is  considered, 
naturally  enabling  numerical  methods  for  computational 
analysis.  The  Ith  element  of  the  k'h  symbol,  spectrally 
sampled  at  uniform  intervals  of  1 oo,  is  given  by 


Fk [l]  =  y]  Af 1 0)o)  e  j°k (Io>° )  =  A[l]  e  J(S>k [,J  (14) 

The  n,h  element  of  the  it  symbol,  temporally  sampled  at 
uniform  intervals  of  t0  such  that  to(Oo  =  2n/Ns  where  A  is  the 
number  of  samples,  is  given  by 


1  Ns-1~ 

fk[n]=—  I  Fk[l]e 

Ns  1  =  0 


j2nln/Ns 


1  N ~  1A[lJej^k[l]+2Kln/Ns 


(15) 


=  —  I 
Ns  1  =  0 


The  sampled  signals  are  further  windowed  in  time  to  be 
realizable.  Windowed  symbol  must  contain  at  least  99.9% 
of  the  original  energy. 
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IV.  Basis  Selection  and  Performance 

To  validate  the  symbol  design  theory  and  the  associated 
metrics,  two  types  of  basis  functions  are  used  to  generate  a 
family  of  M  =  4  symbols  with  intermediate  frequency  (IF) 
bandwidth  of  Qo=  2n  x  103  rad/s  (1.0  KHz)  and  dispersive 
gain  Gd=  10. 

A.  Sinusoidal  Basis 

The  family  of  sinusoidal  functions  is  one  obvious  choice 
of  basis  functions  due  to  their  orthogonal  properties.  Using 
this  family,  the  general  form  for  each  element  in  a  set  of 
sinusoidal  basis  functions  is  described  as: 

( pn(co )  =  y[~2  sin{an  co)  (16) 


These  basis  functions  are  constructed  to  maintain 
orthonormality  within  a  uniform  taper  rather  than  the 
optimal  cosine  taper;  A(co)  in  (9)  is  replaced  with  a  uniform 
taper.  However,  the  symbols  are  still  synthesized  with  the 
cosine  taper  although  the  phase  functions  are  designed  for  a 
uniform  envelope.  This  “mismatch”  is  necessary  to  retain 
the  desirable  sidelobes  characteristics  of  the  cosine  taper 
envelope.  The  ACF  improvement  comes  at  the  cost  of 
possible  degradation  in  the  rms  metric  of  the  CCF.  The 
CCF  metric  will  no  longer  adhere  to  (13)  although  (12)  is 
still  valid  for  the  ACF.  The  correlation  metrics  are  shown  in 
Table  I  while  the  ACF  and  CCF  results  are  shown  Fig  1. 
The  time  function  windowed  at  99.9%  of  the  energy  is  show 
in  Fig  2.  A  comparison  of  Linear  Frequency  Modulated 
(LFM)  symbols  with  equivalent  time  bandwidth  product  is 
also  included  to  provide  a  measure  of  comparison. 


where  ak  e  Integers,  ak  ^  a,  V  k,  l  and  bandwidth  Q0  =  2mn 
for  m ,  a  positive  integer.  These  constraints  ensure 
compliance  to  the  orthonormal  requirement  in  (9).  The 
symbols  were  constructed  using  the  process  outlined  in  (5) 
to  (10).  The  correlation  metrics  of  the  synthesized  symbols 
are  summarized  in  Table  I,  including  Peak  Sidelobe  Level 
(PSL),  Integrated  Sidelobe  Level  (ISL)  and  Peak  Cross- 
Correlation  Level  (PCCL).  The  rms  metrics  of  the  symbols 
meet  the  theoretical  expectation  of  (12)  and  (13)  thus 
validating  the  symbol  design  theory.  The  autocorrelation 
sidelobe  metrics,  PSL  and  ISL,  are  characteristic  of  the 
cosine  taper  envelope  and  remain  relatively  constant  with 
increased  bandwidth,  dispersive  gain,  and  choice  of  other 
basis  functions.  However,  the  PCCL  metric  tends  to  change 
with  dispersive  gain,  bandwidth  and  basis  function 
selection.  As  dispersive  gain  increases,  the  PCCL  metric 
improves  and  the  signal  time  duration  (to  capture  99.9% 
energy)  increases.  Applications  requiring  shorter  duration 
signals  may  not  be  able  to  use  symbols  generated  with  a 
large  dispersive  gain  and  benefit  from  lower  PCCL  levels. 
Thus,  piecewise  linear  basis  functions  are  introduced  to 
“flatten”  CCF  sidelobes  without  the  need  to  increase  the 
dispersive  gain  to  reduce  PCCL’s. 


B.  Piecewise  Linear  Basis 

A  set  of  orthonormal  piecewise  linear  basis  functions  is 
generated  with  each  element  in  this  set  defined  as: 


<Pn(®)  =  Pt 


((O  +  Qq/2)  red 


f  co  +  Q0/2-Qn/8^ 
nn/4 


+  X  (-1 )'"( CO  +  (Om  )  red 
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( „  ..  X 

ft)  -  com 
~QJ2 


+  (co  -Q0/2)  red 


^  co  —  Lin  /  2  +  £2n  /  8  ^ 


*0J_ _ 
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A  =  .  at  =  °l  .  «>„  = 

Q0  cik 


Q„(ii/ak-l) 
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(17) 


V.  NLAS  Performance  Demonstration 

For  concept  demonstration,  four  uniquely  coded,  equal 
energy  symbols  were  generated  using  the  above  piecewise 
basis.  An  Unambiguous  input  signal  sL,(t)  was  created  using 
symbol//!)  and  random  noise  n(t),  as  shown  in  (18),  with 
powers  adjusted  to  achieve  a  -20  dB  signal-to-noise  ratio 
(SNR).  An  Ambiguous  input  signal  sjt)  was  created  by 
adding  two  symbols  to  sjt)  as  shown  in  (19).  The  signal 
powers  in//!)  and //!)  were  adjusted  to  achieve  a  0  dB  SNR 
-  the  ambiguous  signal  powers  are  +20  dB  above  the 
unambiguous  signal  power.  The  remaining  symbol  //!)  was 
used  by  the  NLAS  processor  to  generate  the  Adaptive 
Reserved  Code  Thresholds  (ARCT)  used  for  suppressing 
ambiguous  signal  responses  [3]. 


Sjt)  =fj  t)  +  n(t) 

(18) 

sjt)  =  sjt)  +f2(  t)  +/( t) 

(19) 

Unambiguous  signal  sL,(t)  is  first  input  into  the  NLAS 
processor.  Although  there  is  no  ambiguous  energy  to 
suppress  in  this  case,  the  NLAS  processor  goes  through  all 
suppression  operations  and  provides  a  “colored”  baseline 
for  performance  comparison.  Ambiguous  signal  sA( t)  is  then 
input  into  the  NLAS  processor  and  produces  a  final  output 
as  shown  in  Fig.  4.  A  representative  focusing  and 
suppression  operation  for  one  of  the  undesired  ambiguous 
input  signals  is  shown  in  Fig.  3.  As  illustrated,  all  data 
points  exceeding  the  ARCT  (dashed  line)  are  zeroed  out 
prior  to  subsequent  processing  [3].  The  final  NLAS  output 
response  of  Fig.  4  is  nearly  identical  to  the  “colored” 
response  obtained  from  the  unambiguous  signal  input  case- 
if  complete  suppression  of  f2(t)  and  //!)  occurs,  the 
suppressed  response  would  be  identical  to  the  “colored” 
response.  Despite  the  presence  of  ambiguous  responses  that 
were  +20  dB  higher  than  the  response  of  interest,  the  NLAS 
processor  output  permits  reliable  signal  detection. 
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VI.  Conclusion 

Sinusoidal  basis  functions  provide  a  means  for  generating 
optimal  mutually  dispersive  symbols,  with  dispersive  gain 
providing  control  over  crosscorrelation  dispersion.  The 
sinusoidal  basis  achieves  theoretical  autocorrelation  and 
crosscorrelation  performance,  in  terms  of  rms  time  duration, 
demonstrating  the  validity  of  the  proposed  theorem.  If  the 
envelope  optimality  constraint  is  relaxed  during  phase 
functions  design,  piecewise  linear  basis  functions  provide 
improved  performance  in  terms  of  cross  correlation 
sidelobes  with  minimal  impact  on  the  rms  metric. 

The  theory  of  designing  mutually  dispersive  symbols  and 
the  resultant  generation  process  using  sinusoidal  bases 
provides  a  reliable  method  for  generating  M  mutually 
dispersive  symbols,  with  bandwidth  Qo  controlling 
autocorrelation  compression  and  dispersive  gain  GD 
controlling  crosscorrelation  dispersion.  Although  the 
correlation  rms  metrics  are  invariant  to  basis  function 
selection,  sidelobe  characteristics  of  the  auto-  and  cross¬ 
correlation  functions  are  dependent  upon  the  choice  of 
envelope,  basis  functions,  bandwidth  and  dispersive  gain. 

Using  the  proposed  symbol  design  process  and  piecewise 
basis  functions,  the  ambiguity  suppression  capability  of  a 
NLAS  process  was  demonstrated.  Despite  the  presence  of 
two  ambiguous  signal  responses  that  were  +20  dB  higher 
than  the  response  of  interest,  the  NLAS  process  yielded  an 
output  that  clearly  provides  for  reliable  signal  detection. 

“The  views  expressed  in  this  article  are  those  of  the  author(s)  and 
do  not  reflect  the  official  policy  or  position  of  the  United  States  Air 
Force,  Department  of  Defense,  or  the  US  Government.  ” 
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